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Abstract 



Neutral triple-gauge-boson couplings ZZZ , ZZ'y and Z77 have been stud- 
ied with the DELPHI detector energies between 183 and 
208 GeV. Limits are derived on these couplings from an analysis of the re- 
actions e + e~— > Zj, using data from the final states 7//, with / = q or 1/, from 
e + e - ^ ZZ, using data from the four-fermion final states qqqq, qqn + fi~ , qqe + e~ , 
qqvv, fi + fi~uu and e + e~z/z/, and from e + e~ — > Zj*, in which the final state 7 is 
off mass-shell, using data from the four-fermion final states qqe + e~ and qq^i + . 
No evidence for the presence of such couplings is observed, in agreement with 
the predictions of the Standard Model. 
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1 Introduction 

One of the important properties of the Standard Model which can be tested at LEP2 
is its non-Abelian character, leading to the prediction of triple-gauge-boson couplings. 
However, while non-zero values of these couplings are predicted for the charged (WW'y, 
WWZ) sector, the SU{2) x U{\) symmetry of the Standard Model predicts the absence 
of such couplings in the neutral sector, namely at the ZZZ, ZZ^ and Z77 vertices. This 
paper describes an investigation of this prediction by DELPHI using LEP2 data taken 
between 1997 and 2000 at energies between 183 and 208 GeV. 

1.1 Phenomenology of the neutral triple-gauge-boson vertex 

Within the Standard Model, production of two neutral gauge bosons in e + e~ collisions 
proceeds at lowest order via the t- or -u-channel exchange of an electron. These processes 
are shown in figures la) and b), where both on- and off-shell 7 production is implied, as 
is the subsequent decay of the final state Z or off-shell 7 into a fermion-antifermion pair. 
Figure lc) shows a contribution to production of the same final states which could come 
from physics beyond the Standard Model by the s-channel exchange of a virtual 7 or Z 
via a neutral triple-gauge-boson coupling. In the reactions e + e - ^ Z7 and e + e~^ ZZ 
the final state can, to a good approximation, be considered to be of two on-shell bosons, 
so that only the exchanged boson at the triple-gauge-boson vertex need necessarily be 
considered as off-shell, while in the reaction e + e~^ Z7* both the exchanged boson and 
the outgoing 7* are off-shell. 1 A further process containing a neutral triple-gauge-boson 
coupling with two of the bosons off-shell is shown in figure Id); here a single Z is produced 
in the final state Ze + e~ via fusion of two exchanged vector bosons. 

The phenomenology of the case where two of the three neutral gauge bosons inter- 
acting at the V^V^V® vertex are on mass-shell has been described in [1]. In this case, 
there are twelve independent anomalous couplings satisfying Lorentz invariance and Bose 
symmetry. Calling V the exchanged boson (V = Z,j), the couplings (i =4,5) produce 
a ZZ final state and hf (i — 1 • • -4) the Z7 final state. The couplings h\ and h\ are 
CP-conserving and fY, h\ and h\ are CP- violating. There are no couplings common to 
production of both the ZZ and Z7 final states. 

A complete phenomenological description of the anomalous neutral gauge couplings in 
the case where one, two or three of the gauge bosons interacting at the V^V^V^ vertex 
may be off mass-shell has been developed in [2] . Following the treatment of the charged 
triple-gauge-boson vertex developed, for instance, in [3,4], all the Lorentz-invariant forms 
which can contribute to the ZZZ , ZZj and Z77 vertices are listed, imposing Bose 
symmetry as appropriate. An effective Lagrangian model is then developed in terms of 
the operators of lowest dimension which are required to reconstruct fully all the vertex 
forms, and which affect only the neutral triple-gauge-boson vertex. 2 This leads to a 
Lagrangian with operators of dimension, d, ranging from d = 6 to d = 12. Such an 
expansion is valid in the case where the new physics energy scale, A, represented by the 
operators is very high, at least satisfying the condition A (mz,\/s), and the relative 
contribution from an operator of dimension d may be expected to be suppressed by 
a factor l/A^ d ~ 4 \ In the analysis we report here, we have considered only the lowest 

throughout this paper, we write "V*" when we wish to be explicit that a vector boson V is off mass-shell. When it is 
clear that it is on mass-shell, or when it can be either on or off mass-shell, the star ("*") is omitted. 

2 The V^VjVg vertex functions receive contributions from both transverse and scalar terms, the latter contributing in 
the case where one off-shell Z decays to a heavy fermion pair through its axial coupling. In the analysis of LEP data only 
transverse terms need be considered, due to the negligible contribution of Z —*tt decays. The contribution of scalar terms 
is therefore ignored in the following. 



a) 




Figure 1: a), b) Lowest order Feynman diagrams for the production of two gauge bosons 
ZZ and Zq by Standard Model processes, where both on- and off-shell 7 production is 
implied, c) Production of the same final states via an anomalous interaction among three 
neutral gauge bosons, d) Production of the Ze + e~ final state via an anomalous neutral 
triple-gauge-boson coupling. 



dimension operators contributing to the parameters we have determined. In addition to 
satisfying Lorentz and Bose symmetry, the operators are required to be 6 r (l) em -invariant, 
and both CP-conserving operators, O, and CP-violating operators, O, with coefficients 
£ and £, respectively, are considered: 

c = e(J2 tT° v§ oT° v§ + E ~^dT M ) ■ (i) 

i,CP+ i,CP- 

Of the operators included in the sum defined above, some affect the V°ZZ and V Z , y* 
vertices (V° = 7), some the V°Z / y* vertex only, and some the V°Z / y* and V°Z r y 
vertices; none contribute to all three vertices. In [2] a connection is made between the 
coefficients l{ and li of the operators in the effective Lagrangian describing the general 
\1V2V3 vertex and the dimensionless coefficients h and / describing on-shell Zy and ZZ 
production, respectively: retaining only the terms corresponding to contributions from 
operators of lowest dimension, each of the coefficients h\ , h% , fX and fY (which are 
dimensionless) is related to one operator of dimension d = 6 by /, h = £ Vl v z y s m%. The 
lowest dimension operators contributing to h\ and h% have d = 8. 

As in the case of the charged triple-gauge-boson couplings, a further simplification in 
the possible structure of the effective Lagrangian is obtained by imposition of SU{2) x 
£7(1) invariance on its form. Such a form is presented in an addendum [5] to [2], and 
the effective Lagrangian reduces to a sum of two terms, both with dimension d — 8, 
one (Osu(2) x u(i)) CP-conserving and one (Osu(2)xU(i)) CP-violating. This simplification 



leads to constraints between some of the 



v°v°v° 



or 



n VPVnV? 



defined in equation (1) 
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-j4r/(2)x[/(i) , (2) 
v 2 ~ 

-j£sU(2)xU(l) , (3) 

where $w is the weak_ mixing angle, v is the vacuum expectation value of the Higgs 
field and £su(2)xu(i), £su(2)xu(i) are the coefficients of the operators Osu(2)xU(i) and 
Osu(2)xU{i), respectively. If applied solely to the on-shell channels Z7 and ZZ, these 
conditions become, respectively: 

v 2 

f z cotO w = h z = -jl = -hjtan6 w = m 2 z —£ su{2)xU{1) , (4) 

v 2 ~ 

f z cot 6 W = h z = -f] = -hj tan 9 W = m 2 z —£ su(2)xU(1) . (5) 

The SU(2) x £/(I)-conserving Lagrangian considered in [2] is constructed so as to affect 
only the neutral gauge boson and Higgs sectors, and an alternative form, which would 
additionally affect off-mass-shell charged gauge boson production, has been proposed 
in [6]. This leads to a Lagrangian with four possible terms, two CP-conserving (0$, 0%) 
and two CP-violating (Of, Og) and hence to looser constraints between the possible 
contributing operators: in each of the sets of conditions (2) - (5) listed above, the relations 
corresponding in the diagrams of figure 1 to 7 and Z exchange decouple, giving, for 
instance in the case of (2), the separate conditions 

£ zzz cot 6 W = £ ZZl and £ 2 ZZ7 = £^ tan 6 W , (6) 

with an analogous separation in conditions (3) - (5). This leads to four coefficients, £g' B , 
£%' B , related to the respective Lagrangian operators by appropriate factors of mz and 
v. In both the stronger and weaker of these sets of constraints (which we refer to as 
the Gounaris-Layssac-Renard and Alcaraz constraints, with respect to the authorship 
of references [2] and [6]), the gauge- invariant operators all now contribute to all three 
neutral triple-gauge-boson vertices, V°ZZ, V°Z^* and V°Z^. 

In order to study the V^V^Vg vertex, three physical final states have been defined 
from the data: Zj, Z7* and ZZ. The first of these is a three-body final state comprising 
the Z decay products and a detected photon, while the other two are four-fermion final 
states with, respectively, one or two fermion-antifermion pairs having mass in the Z 
region. Given the phenomenology summarized above, we have then chosen to determine 
the following parameters in our study: 

• Using data from the final states ZZ and Zrf*, values are determined for the coefficients 
of each of the four d = 6 operators which are related in the on-shell limit to one 
of the / coefficients defined in the on-shell formalism of reference [I]. Similarly, 
using data from the final states Zy, Z7* and ZZ, values of the coefficients of the four 
d = 6 operators related to the on-shell h coefficients are determined. In the latter 
case, the ZZ data are used as well as the Z7*, as the off-shell 7 couples to the // 
system over the whole of the four-fermion phase space. However, in these studies the 
statistical contribution of the off-shell final states compared to that of on-shell Z7 or 
ZZ production is very small, so that the values of the parameters determined, quoted 
in dimensionless form, v z v s m 2 z , are directly comparable with published results 
using data from on-shell channels, and the relevant respective likelihood distributions 
may be combined. 



e^cotOw = £ x ZZl = -£ 2 ZZ ^ = -^'tanO 



w 



\ zzz cot9 w = = -£ ZZl = -4 Z77 tan0 



w 
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• The V Z7* vertex is studied on its own by determining the coefficients of the lowest 
dimension operators which affect solely these vertices. There are two such operators, 
both of dimension d — 8, one CP-conserving and involving s-channel 7 exchange 
in the production process illustrated in figure lc), and the other CP- violating and 
involving s-channel Z exchange in the same diagram. Again, data from both the 
Zy* and ZZ final states were used in the determination of the coefficients of these 
operators, and the coefficients are quoted in dimensionless form: t v ° v ° v °m%. 

• The coefficients of the SU{2) x £/(l)-conserving operators are determined, using 
both the Gounaris-Layssac-Renard and the Alcaraz constraints. They are quoted in 
a dimensionless form, such that in the on-shell limit they become equal to one of 
the hf occurring in the constraint equations (4) and (5) above. 

A list of the parameters we have determined, the definitions of the operators to which 
they refer and (where relevant) the on-shell coefficients to which they are related is given 
in table 1. 

1.2 Experimental considerations 

Of the three final state channels, Zy, ZZ and Z7*, defined in the previous section, 
the most precise limits on anomalous couplings are derived from the first, when the 
final state photon is on-shell. In this channel, the kinematic region with high photon 
energy and large photon polar angle is most sensitive to the anomalous couplings, and in 
this region the anomalous interactions give rise to a change in the total rate and to an 
enhancement of the production of longitudinally polarized Z bosons. Our analysis covers 
two reactions to which the diagrams describing Zy production provide the dominant 
contribution: e + e~ — > z/z/7, in which the observed number of events is compared with 
the number predicted from the total cross-section for this process, and e + e~ — > qqy, with 
the qq system coming predominantly from Z decay, in which the distribution of the decay 
angle of the Z in its rest frame with respect to the direction of the Z in the overall centre 
of mass is studied. The present analysis uses data from LEP2 at energies ranging from 
189 to 208 GeV. Previous DELPHI results on this channel can be found in [7]; they used 
data with energies up to y/s = 172 GeV, and the limits were obtained using an analysis 
based only on the value of the observed total cross-section. 

The total ZZ cross-section is also sensitive to the anomalous couplings, and the sen- 
sitivity increases strongly with y/s. Large interference between Standard Model and 
anomalous amplitudes arises for CP-conserving couplings (especially for /rf) when con- 
sidering the differential cross-section da/ d | cos 6z\, where 9 Z is the Z production angle 
with respect to the beam axis. The analysis reported here is based on a study of this 
differential distribution in the LEP2 data in the energy range 183 to 208 GeV. DELPHI 
has previously reported a study of the ZZ - production cross-section in all visible /1/2/3/4 
final states in these data [8]. The same sets of identified events have been used in the 
present analysis, with the exception of the qqr + r~ , t + t~vu and l + l~l + l~ channels, which 
are not used. 

In a separate publication [9], DELPHI has studied the Zy* final state in the same 
LEP2 data as used for the channels described above, reporting on a comparison of the 
cross-section for Zy* production in various channels with Standard Model predictions. 
We use the samples identified in [9] in the qqe + e~ and qqfi + fi~ final states in the present 
analysis, which thus represents an interpretation of these data for the first time in terms 
of possible anomalous gauge couplings. The data were examined as a function of the 
bidimensional (iWj+j-, M q q) mass distribution, requiring one of them to be in the region 
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Table 1: Parameters determined in this study, corresponding Lagrangian operators in 
the models of references [2] and [6] , and (where appropriate) related on-shell parameters: 
a) Coefficients of lowest dimension operators contributing to ZZ and Zy* production or to 
Zy* and Zy production; b) Coefficients of lowest dimension operators affecting only the 
V°Z'~f* vertices; c) Coefficients of SU(2) x £/(l)-conserving operators according to i) the 
Gounaris-Layssac-Renard constraints and ii) the Alcaraz constraints. The constraints are 
given in the text. The vertices V^V^V^ affected by these operators (without distinguish- 
ing the V t ° as on- or off-mass-shell) are indicated in column 1. The fields Z M , F^ B p and 
Wfj, represent the Z, photon, U(l)y and SU(2) L fields, respectively; Z MI/ , F pu and B pu are 
the contractions of the respective field tensors with the four-dimensional antisymmetric 
tensor; <£> is the Higgs field and v its vacuum expectation value, D represents the covariant 
derivative of SU(2) x £7(1), and 17 are the Pauli matrices. 
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of the Z mass, and they were also divided into two regions of the l + l~ polar angle with 
respect to the beam direction (equivalent to the variable 6 Z used in the analysis of ZZ 
events). 

Limits on anomalous neutral gauge couplings in the Z7 and ZZ final states have been 
reported by other LEP experiments; recent published results may be found in the papers 
listed in [10,11]. 

2 Experimental details and analysis 

Events were recorded in the DELPHI detector. Detailed descriptions of the DELPHI 
components can be found in [12] and the description of its performance and of the lumi- 
nosity monitor can be found in [13]. The trigger system is described in [14]. For LEP2 
operations, the vertex detector was upgraded [15], and a set of scintillation counters was 
added to veto photons in the blind regions of the electromagnetic calorimetry at polar an- 
gles around 9 = 40°. The performance of the detector was simulated using the program 
DELSIM [13], which was interfaced to the programs used in the generation of Monte 
Carlo events and to the programs used to simulate the hadronization of quarks from Z 
and 7* decay or from background processes. During the year 2000, one sector (1/12) of 
the Time Projection Chamber (TPC), DELPHI'S main tracking device, was inactive for 
about a quarter of the data-taking period. The effect of this was taken into account in 
the detector simulation and in the determination of cross-sections from the data. 

The selection of events in the three physical final states, Z7, ZZ and Z^y*, considered in 
this paper, and the simulation of the processes contributing to signals and backgrounds 
are described in the following subsections. In the case of the ZZ and Z7* samples, the 
reader is referred to recent DELPHI publications on the production of these final states 
(references [8,9], respectively) for a full description of the event selection procedures. The 
event samples used in the present analysis of these two final states have been selected 
using essentially identical procedures to those described in [8,9], and cover the same energy 
range (183 - 208 GeV). These procedures are summarized, respectively, in sections 2.2 
and 2.3 below, and any changes from the methods described in [8,9] are mentioned. 
DELPHI has also reported a study of events observed at LEP2 in which only photons 
and missing energy were detected [16]. The present analysis uses data in the part of the 
kinematic region covered in [16] in which a high energy photon is produced at a large 
angle with respect to the beam direction; data in the energy range 189 - 208 GeV have 
been used. The selection procedures specific to this final state as well as to that in which a 
quark- ant iquark pair is produced, rather than missing energy, are described in section 2.1 
below. 

In the final year of LEP running, data were taken over a range of energies from 205 to 
208 GeV. The values of the centre-of-mass energy quoted in the descriptions below for 
that year correspond to the averages for the data samples collected. 

2.1 The Z7 final state 

The selection procedure for Zq production in the kinematic region with greatest sen- 
sitivity to anomalous gauge couplings concentrated on a search for events with a very 
energetic photon in the angular range 45° < # 7 < 135°, where 6> 7 is the polar angle of the 
photon with respect to the beam direction. This angular region is covered by DELPHI'S 
barrel electromagnetic calorimeter, the High density Projection Chamber (HPC). The 
search was conducted in events with two final state topologies: vv^ and qcpf. 
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The vv^ sample was selected from events with a detected final state containing only a 
single photon. Its energy, E 7 , was required to be greater than 50 GeV and only photons 
in the range covered by the HPC, 45° < 9 y < 135° were accepted. No tracks or hits were 
allowed in the TPC. It was also required that no electromagnetic showers with energy 
exceeding defined background noise levels were present in the forward electromagnetic 
calorimeter and the luminosity monitor. Further showers in the HPC were accepted only 
if they were within 20° of the first one, and such showers were then combined. Cosmic 
ray events were suppressed by requiring that any signal in the hadronic calorimeter be 
in the same angular region as the signal in the electromagnetic calorimeter, and that the 
electromagnetic shower point towards the beam collision point within an angle of 15°. 
The trigger efficiency was measured using Compton and Bhabha events. The expected 
numbers of events were calculated using the generators NUNUGPV, based on [17], and 
KORALZ [18], interfaced to the full DELPHI simulation program. The results obtained 
applying these criteria are shown in table 2. From the simulations, the efficiency for 
detecting z/z/7 events in the kinematic region considered here was shown to be independent 
of the centre-of-mass energy, with an average value of (50.7 ± 2.0)% for the data sample 
listed in table 2. Contributions from background sources to this channel were estimated 
to be negligible. Full details of the analysis of this final state may be found in [16]. 
The distribution of x 1 , the energy of identified photons normalized to the beam energy 
(x 7 = E^/Ebeam) before imposing the cut at E 1 = 50 GeV is shown for photons with 
x 7 > 0.05 in figure 2a), which also shows the expectation of the Standard Model. 





Integrated 


Selected 


Total predicted 


(GeV) 


luminosity (pb _1 ) 


data 


events 


188.6 


154.7 


87 


89.2 


191.6 


25.1 


14 


13.1 


195.5 


76.2 


32 


37.5 


199.5 


83.1 


45 


38.5 


201.6 


40.6 


20 


18.2 


206.1 


214.6 


98 


102.3 


Total 


594.3 


296 


298.8 



Table 2: z/z/7 fi na l state: Integrated luminosity and numbers of observed and expected 
events at each energy, s/s. 

In the selection of events in the q<pi channel, the same requirements were imposed 
on the energy and polar angle of photon candidates as in the z/z/7 case, namely: E 1 > 
50 GeV and 45° < # 7 < 135°. In addition, events were required to have six or more 
charged particle tracks, each with length greater than 20 cm, momentum greater than 
200 MeV/c, polar angle between 10° and 170°, and transverse and longitudinal impact 
parameters at the interaction point of less than 4 cm. The total charged energy in 
the event was required to exceed 0.10^ and the effective energy of the collision [19], 
excluding the detected photon, v^s 7 , was required to satisfy \fs i < 130 GeV. Jets were 
reconstructed using the LUCLUS [20] algorithm and, omitting the 7, the event was forced 
into a two-jet configuration. The identified photon was required to be isolated from the 
nearest jet axis by at least 20°. The efficiency, purity and the expected numbers of 
events from qq(^j) production were computed using events generated with PYTHIA [20], 
relying on JETSET 7.4 [20] for quark fragmentation, and interfaced to the full DELPHI 
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simulation program. The results obtained applying this procedure are shown in table 3. 
The efficiency for detecting qq'j events in the kinematic region considered here was found 
to be almost independent of the centre-of-mass energy for the data sample used, with an 
average value of (76.4 ± 0.2)%. The main background processes, contributing about 3% 
of the sample, came from qq production with a photon from fragmentation of one of the 
quarks, and from WW production. 



(GeV) 


Integrated 
luminosity (pb _1 ) 


Selected 
data 


Total predicted 
events 


Expected 
background 


188.6 


154.3 


454 


467.3 


14.9 


191.6 


25.4 


79 


75.0 


2.6 


195.5 


77.1 


203 


214.1 


5.8 


199.5 


84.2 


208 


225.5 


5.9 


201.6 


40.6 


130 


104.5 


2.8 


205.9 


218.8 


507 


515.1 


13.9 


Total 


600.4 


1581 


1601.5 


45.9 



Table 3: qq'j final state: Integrated luminosity, numbers of observed and expected events 
and predicted background contribution at each energy, y/s. 



Summing over all energy points, totals of 296 and 1581 events were observed in the 
vV'-f and qq'j channels, respectively. These numbers may be compared with the totals 
expected from simulated production of these final states by Standard Model processes: 
298.8 events in uuj, and 1601.5 events in qqj. 

In the z/X/7 channel, values of the gauge boson coupling parameters were derived by 
comparing the observed number of events with that predicted from the total cross-section 
for this process, while in the qq'j channel a fit was performed to the distribution of | cos a*\, 
where a* is the angle of the quark or antiquark from Z decay in the Z rest frame with 
respect to the direction of the Z in the overall centre of mass. The value of | cosa:*| was 
estimated from the directions of the vectors in the laboratory frame p 7 and pi (i — 1, 2) 
of the reconstructed photon and jets, respectively, from the relation: 

cot a* = 7^ cotai — — V (7) 

V smai/ 

a sin («l + «2) P 7 Pi , 1 fa , 

with p = , cos a; = —-, — — ; — r and 7 = = . (8) 

sm «i + sm a 2 IP7I' |Pi| V 1 ~~ P 

The distribution of | coso;*| for the data selected in the qq'j channel is shown in figure 2b) 
and compared with the predictions of the Standard Model and of a non-standard scenario 
with hj = ±0.2. The predictions for non-zero neutral gauge boson couplings were made 
by reweighting the simulated samples produced according to the Standard Model with 
the calculations of Baur and Berger [21] 3 . 



2.2 The ZZ final state 

The study of the triple-gauge-boson vertex in ZZ production used the samples of 
events selected in the qqqq, qq^ + ^~ , qqe + e~ , qquu, fi + fi~i>v and e + e~z/z/ final states. The 
procedures used to extract the data have been described fully in [8]; we give here a brief 

3 The code used was modified by a factor i according to the correction suggested by Gounaris et al [22], 
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summary of the methods used in the selection of events in each of these final states, and 
provide a table of the total numbers of events observed and expected for production of 
each of them by Standard Model processes. 

The ZZ — > qqqq process represents 49% of the ZZ decay topologies and produces four 
or more jets in the final state. After a four-jet preselection, the ZZ signal was identified 
within the large background from WW and qqj production by evaluating a probability 
that each event came from ZZ production, based on invariant-mass information, on the 
6-tag probability per jet and on topological information. 

The process e + e~ — > qql + l~ has a branching ratio in production of 4.7% per lepton 
flavour. High efficiency and high purity were attained with a cut-based analysis using the 
clear experimental signature provided by the two leptons, which are typically well isolated 
from all other particles. The on-shell ZZ sample was selected by applying simultaneous 
cuts on the masses of the / + /~ pair, on the remaining hadronic system and on their sum. 

The decay mode qqvv represents 28% of the ZZ final states. The signature of this decay 
mode is a pair of jets, acoplanar with respect to the beam axis, with visible and recoil 
masses compatible with the Z mass. The most difficult backgrounds arise from single 
resonant Wev e production, from WW production where one of the W bosons decays into 
tv t , and from qq events accompanied by energetic isolated photons escaping detection. 
The selection of events was made using a combined discriminant variable obtained with 
an Iterative Discriminant Analysis program (IDA) [23]. 

The final state l + l~vv has a branching ratio in ZZ production of 1.3% per charged 
lepton flavour. Events with I = fj,,e were selected with a sequential cut-based analysis. 
The on-shell ZZ sample was selected from the events assigned to this final state by 
applying cuts on the masses of the l + l~ pair and of the system recoiling against it. 
The most significant background in the sample is from WW production with both Ws 
decaying leptonically. 

In the estimation of the expected numbers of events in all the final states discussed 
above, processes leading to a four-fermion final state were simulated with EXCAL- 
IBUR [24], with JETSET 7.4 used for quark fragmentation. Amongst the background 
processes leading to the final-state toplogies described above, GRC4F [25] was used to 
simulate Wev production, PYTHIA for 55(7), KORALZ for n + fi~( / y) and t + 1-^(7), BH- 
WIDE [26] for e + e"(7), and TWOGAM [27] and BDK [28] for two-photon processes. 

The presence of anomalous neutral triple-gauge-boson couplings in the data samples 
described above was investigated by studying the distribution of the Z production polar 
angle, | cos8z\- For events in the qqvv and l + l~vv final states, the Z direction was taken 
to be the direction of the reconstructed di-jet or pair, respectively, while in the 
qql + l~ final state, the Z direction was evaluated following a 4-constraint kinematic fit to 
the jet and lepton momenta, imposing four-momentum conservation. In the qqqq final 
state, the indistinguishability of the jets leads to three possible jet-jet pairs, each of which 
could come from ZZ decay. A 5-constraint kinematic fit was performed on each of these 
combinations, imposing four-momentum conservation and equality of the masses of the 
two jet pairs. The fit with the minimum value of \ 2 was retained and the value of | cos 9z\ 
evaluated from the fitted jet directions. 

Figure 3 shows the distribution of | cos8z\ for a high purity sample of ZZ data, com- 
posed of the qql + l~ and l + l~vv samples defined above and, for illustrative purposes, 
samples of qqqq and qqvv events defined by stringent cuts on the probabilistic variables 
used in these channels (qqqq probability > 0.55, and qqvv IDA variable > 3), so as to 
suppress the background levels present. (As described below, no cuts were imposed on 
these variables in the determination of coupling parameters). The figure also shows the 
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Standard Model expectations and the distributions predicted for values of /if = ±1.5. 
The content of this sample is shown in table 4. The selection efficiencies for all of the 
channels analysed were found to have little energy variation over the range considered 
here, and the values quoted are averages for the whole experimental sample. They repre- 
sent the fraction of events of the relevant four-fermion final state present in the selected 
sample, while the estimated backgrounds refer to contributions from other channels. 



Channel 


Integrated 


Selected 


Total predicted 


Expected 


Selection 




luminosity (pb _1 ) 


data 


events 


background 


efficiency 


qqqq 


665.1 


76 


69.4 


22.1 


0.18 


qqH + H~ 


665.3 


21 


22.0 


1.1 


0.86 


qqe + e~ 


665.3 


19 


23.7 


2.6 


0.73 


qqvv 


639.0 


45 


55.5 


22.3 


0.21 


m-w 


665.3 


10 


8.9 


4.7 


0.30 


Total 




171 


179.5 


52.8 





Table 4: ZZ production: Integrated luminosity, numbers of observed and expected events 
and predicted background contribution for each topological final state, summed over all 
energies. The last column shows the energy-averaged selection efficiency for each final 
state. 

In the determination of the coupling parameters, extended maximum-likelihood fits 
were made to the distribution of | cos 9z\ for data from the channels selected with cut- 
based analyses (qql + l~ and l + l~vv), while for the channels selected using probabilistic 
methods (qqqq and qqvv), a simultaneous fit was made to the distributions of |cos#z| 
and of the discriminant variable (the ZZ probability for qqqq and the IDA output variable 
for qqvv), without applying any cuts on the values of these variables. 

The predictions for non-zero neutral gauge boson couplings were made by reweighting 
the simulated samples produced according to the Standard Model with the calculations of 
the DELTGC [29] event generator, which adds the amplitude from hypothesized neutral 
triple-gauge-boson vertices to all the other amplitudes contributing to the production of 
any four-fermion final state. 

2.3 The Z 7 * final state 

In a separate publication [9], DELPHI has reported on a study of Z7* production in 
LEP2 data, and in particular on a comparison of the observed cross-section with Stan- 
dard Model predictions, using data from a variety of four-fermion final state topologies 
involving both hadronic and leptonic Z decay modes. In the present analysis, we interpret 
data in the qqn + n~ and qqe + e~ final states in terms of possible anomalous triple-gauge- 
boson interactions. These two channels are chosen because the two final state leptons are 
typically well isolated from all other particles, allowing such events to be selected with 
high efficiency over the whole region of 7* mass. Events with either the / + /~ or the qq 
invariant mass in the vicinity of the Z mass and the other invariant mass not in the Z 
region were then used in the estimation of possible anomalous gauge coupling parameters. 
Full details of the selection procedure are given in [9]; a summary of the main features 
follows. 

Events containing total charged hadronic energy above 0.3\/i and at least two lepton 
candidates of the same flavour and opposite charge were selected. All particles except the 
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lepton candidates were clustered into jets and a kinematic fit requiring four-momentum 
conservation was applied. At least one of the two lepton candidates was required to 
satisfy strong lepton identification criteria, while softer requirements were specified for 
the second. In order to increase the purity of the selection, further cuts were made in 
two discriminating variables: P t mm ; the lesser of the transverse momenta of the lepton 
candidates with respect to their nearest jet, and the x 2 P er degree of freedom of the 
kinematic fit. This procedure selected a total of 170 events in the combined qqn + n~ and 
qqe + e~ channels. The Z7* sample was then defined within the selected qql + l~ data by 
requiring the mass of one and only one // pair to be in the Z region. This was effected by 
imposing mass cuts in the (M hadrons , M a+n -) and (M hadrons , M e+e -) planes, where M hadrons 
represents the mass of the qq pair estimated from the reconstructed hadronic data; these 
cuts are defined in figures 4a) and b) for the qqn + f^~ and qqe + e~ samples, respectively. 
Table 5 summarizes the selection procedures outlined above, showing, for the sum of 
data over all energy points, the total integrated luminosity, the numbers of observed 
and predicted events in the Z7* region, defined as described above, and the estimated 
selection efficiency (defined as for the ZZ sample described in section 2.2 above) for each 
topological final state. The selection efficiencies for the qqn + fi~ and qqe + e~ samples 
analyzed here were shown to be independent of the centre-of-mass energy, with average 
values of (44.1 ± 0.2)% and (28.7 ± 0.1)%, respectively. The backgrounds in the selected 
samples are small, coming mainly from qqr + r~, WW and, in the case of qqe + e~, from 
qq{l) production. In the estimation of backgrounds and selection efficiency, the simulation 
of processes leading to four-fermion final states was done with WPHACT [30] , using the 
JETSET model for quark hadronization, while the qql'j) final state was simulated with 
the KK2f [31] model. Both of these programs were interfaced to the DELPHI simulation 
program. 



Channel 


Integrated 
luminosity (pb _1 ) 


Selected 
data 


Total predicted 
events 


Expected 
background 


Selection 
efficiency 


qqfi + ^- 


666.7 


35 


36.7 


3.4 


0.44 


qqe + e~ 


666.7 


39 


36.3 


6.0 


0.29 



Table 5: Zy* production: Integrated luminosity, numbers of observed and expected events 
and predicted background contribution for each topological final state, summed over all 
energies. The last column shows the energy-averaged selection efficiency for each final 
state. 

Some aspects of the phenomenology of qql + l~ production in the context of possible 
neutral triple-gauge-boson couplings, and of the data selected in the qqfi + fi>~ and qqe + e~ 
channels, are demonstrated in figures 4 and 5. Figures 4a) and b) show the distribu- 
tions in the (Mhadrons, M t +i~) planes predicted by the Standard Model for the qqfi + fi~ 
and qqe + e~ final states, respectively. These differ considerably, due to the presence of 
additional diagrams contributing to qqe + e~ production, in particular those correspond- 
ing to the production of Ze + e~ and 7*e + e _ by t-channel processes. These effects have 
been discussed fully in [9]. The effect of an anomalous triple-gauge-boson coupling in 
these channels is illustrated in figures 4c) and d), which show, respectively, the difference 
between the expected distributions on the (M qq: M^+a-) and (M qq , M e + e -) planes when 
a non-zero contribution from the d = 8 operator 4 ZZl (defined in table 1) is included, 
and when only the Standard Model amplitudes are used. Again, some differences be- 
tween the predictions for the qqii + \x~ and qqe + e~ final states are observed; these are due 
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to the presence of additional diagrams in the qqe + e~ amplitude, in this case the V V 
fusion diagram leading to Ze + e~ production, shown in figure Id). The overall effect is a 
negative interference between s- and t-channel amplitudes: for the example shown, the 
predicted content of figure 4d) (qqe + e~) is ~ 40% of that of the qqfi + fi~ prediction. 

Data selected over the whole region of the qql + l~ phase space are presented in fig- 
ures 5a) and b) in the form of the distributions of Mi+i- (I = fi, e) and M qq . These plots 
also show the expectations of the Standard Model and of a model in which an anomalous 
contribution l/^m^ = 3-4 from the operator 4 ZZl is present. 

In the determination of the coupling parameters, the regions in the plane of the masses 
of the two fermion-antifermion pairs defining the Z^y* samples in the qq^ + ^~ and qqe + e~ 
final states, shown in figures 4a) and b), respectively, were divided into a small number 
of bins of unequal size, but containing roughly equal numbers of events predicted by the 
Standard Model. Different bin definitions were made for the two channels; the bins are 
also defined in the figures, and they correspond to those used by DELPHI in [9] in the 
determination of the Z7* cross-section. In [9], each of the mass bins defined for the qqe + e~ 
event sample was further divided into two angular regions, (40° < 9i+ t - < 140°) and 
(9i+i- < 40° or 81+1- > 140°), where &i+i- is the polar angle of the final state e + e~ system 
with respect to the beam direction. These angular regions correspond to DELPHI'S barrel 
and endcap regions, respectively. In the present analysis, we have extended this division 
to apply to muon as well as electron pairs in the qql + l~ final states. Binned likelihood fits 
to the couplings were then made with the bins in (M qq , M1+1-) and 9i+ t - thus defined. As 
in the case of the ZZ final state previously described, the predictions for non-zero neutral 
gauge boson couplings in the Z7* data were made by reweighting the simulated samples 
produced according to the Standard Model with the calculations of DELTGC. 

3 Results 

In this section the results of our study are presented, expressed in terms of the param- 
eters listed in table 1 describing the neutral triple-gauge-boson effective Lagrangian. In 
summary, these parameters represent: 

a) the coefficients of the lowest dimension operators contributing to production either 
of the ZZ and Z7* final states, or to production of the Zj, Z^y* and ZZ final states; 
in the on-shell ZZ or Z7 limit each of these parameters becomes equal to one of the 
on-shell coefficients fY or hf; 

b) the coefficients of the lowest dimension operators affecting only the V°Z / f* vertex; 

c) the coefficients of the £77(2) x £/(l)-conserving operators describing the V^V^V^ 
vertex in i) the Gounaris-Layssac-Renard and ii) the Alcaraz formulations. 

The labels a), b), c) above refer to table 1. 

Limits on the parameters at the 95% confidence level are given in table 6 and the 
corresponding likelihood curves are shown in figures 6-10. In all cases, the values quoted 
are derived from one-parameter fits to the data in the Zj, ZZ and Z7* channels described 
in sections 2.1, 2.2 and 2.3 above, summing the distributions from different channels 
where appropriate. In each fit, the values of the other parameters were set to zero, their 
Standard Model value. The results shown include contributions from both statistical and 
systematic effects. 

For reference, we summarize here the composition of the likelihood function from each 
of the final states used in the analysis, described in more detail in the sections above: 
In the Z'-f — > 1/1/7 channel, the number of events with a high energy photon emitted at 
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large polar angle was used in the fit, while in the Z7— > qcpi channel the fit was performed 
to the distribution of the decay angle of the Z in its rest frame. In the channels ZZ — > 
qql + l~ and ZZ — > l + l~vu the distribution of the Z production angle was fitted; in ZZ 
— > qqqq and ZZ — > ggz/z/ simultaneous fits were made to the Z production angle and, 
respectively, to the event probability or discriminant variable distributions. In the Z7* 
channels studied (Z7* — > qqfi + fi~ and Z7* — > qqe + e~) the likelihood was evaluated in 
bins of qq or / + /~ mass and of the polar angle of the detected l + l~ system. 

It may be noted that, in the models conventionally used to describe anomalous gauge- 
boson couplings, including the one used in this paper, all observables have a quadratic 
dependence on the fitted parameters. This effect, which has been previously noted (see, 
for example, [10]), leads to log-likelihood distributions which can have double mimima, 
asymmetries, and a broadening compared with that expected in the Gaussian case. Such 
features are indeed seen in several of the plots in figures 6-10. The confidence limits 
reported in table 6 must therefore be interpreted with this effect in mind. 

3.1 Systematic errors 

In the determination of the confidence limits shown in table 6 and the likelihood curves 
of figures 6-10, several sources of systematic error were considered for each of the final 
states included in the study. These are described below. 

In the vV^f and qcpi channels contributing to Z7 production, uncertainties of ±1% were 
assumed in the values assumed for the Standard Model production cross-sections [17,18], 
and an experimental uncertainty of ±1% was assumed for the energy calibration of the 
electromagnetic calorimeter. The effect of an uncertainty of ±1% in the luminosity mea- 
surement was also computed, while the uncertainties in the calculations arising from the 
finite simulated statistics in signal and background channels and from the uncertainty 
in the knowledge of the background cross-section were found to be negligible in both 
channels. In the vv^j channel, the error due to the uncertainty of ±3% in the trigger 
efficiency was included. In the qcpi channel, the uncertainty in the use of PYTHIA as 
the hadronization model was taken into account by comparing events simulated with 
PYTHIA and HERWIG [32] ; this gave rise to an estimated systematic error on the selec- 
tion efficiency of ±1.7% from this source. In the combination of data at different energies, 
all the above effects were considered as correlated. The resulting overall systematic error 
in the coupling parameters was found to be of the order of 30% of the statistical errors in 
the case of hf and /if, about 50% of the statistical error for K[, and of the same order as 
the statistical error for hj. In combination with Z7* and ZZ data to produce the limits 
on the parameters £ 1 ZZ7 m|, £ l zz "'m 2 z , l Zll mz and £ 1 Z77 to| shown in table 6a) the Z7 
data dominate (see sections 1.1 and 3.2 for further discussion of this point), so that the 
ratios of systematic to statistical errors quoted above are also applicable to the respective 

ij 1 V2 Vs m 2 z results reported in the table. 

A full description of the treatment of systematic effects in the channels contributing 
to production has been given in [8]. In the qqqq channel, the dominant effect arises 
from uncertainties in the modelling of the main source of background, namely production 
of the qq{^) final state, when the subsequent hadronization of the quarks gives rise to 
several jets. In the present study, the effect of this background was estimated by assuming 
an uncertainty of ±5% in the 99(7) production cross-section. In the qql + l~ channel, the 
dominant systematic effect relevant to the present study comes from the uncertainty in 
the efficiency for selecting qqe + e~ and qqfj, + \i~ events, taken to be ±3%. In addition, 
in the qqe + e~ channel a relative uncertainty of ±15% was estimated in the calculation 
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Parameter 


Channels 


95% Confidence 


Related on-shell 




used 


interval 


coefficient 


a) 


£^ Z7 "m 2 z 

£ ZZ V Z 
£ 2 ZZ V Z 


ZZ Zy* 
ZZZY 

zzzy 

ZZZY 


[-0.40, +0.42] 
[-0.38, +0.62] 
[-0.23, +0.25] 
[-0.52, +0.48] 


fi 
fi 
fi 
fi 


l{' Zl m\ 

^ zz v z 

£ z ^m 2 z 


Zq ZY ZZ 
Z 7 ZY ZZ 

z 7 zy zz 


[-0.23, +0.23] 
[-0.30, +0.16] 
[-0.14, +0.14] 


hf 
h! 
hi 




Zj zy zz 


[-0.049, +0.044] 


hi 


b) 


t2 Zrr < 


zzzy 

ZZZY 


[-1.67, +1.92] 
[-0.49, +0.61] 




c) i) 


— cot 0wTn z ^£ S u(2)xu(i) 

- cot w m 2 z ^-£ S u(2)xU(i) 


Zq ZY ZZ 


[-0.13, +0.13] 


hi 


z 7 zy zz 


[-0.045, +0.047] 


hi 


ii) 


-cot 9 w m 2 z ^£f 
-cot 9 w m 2 z ^£$ 


z 7 zy zz 


[-0.14, +0.14] 


hi 


z 7 zy zz 


[-0.049, +0.045] 


hi 


- cot 9 w m%Yf 
-cot O w m 2 z \l$ 


z 7 zy zz 


[-0.23, +0.24] 


hf 


z 7 z 7 * 


[-0.30, +0.18] 


hi 



Table 6: Results of the study of neutral gauge couplings. For each of the parameters 
listed in table 1, the table shows the experimental channels used and the 95% confidence 
limits obtained. The right-most column indicates the parameter which, in the on-shell 
limit, is equal to the parameter determined. In the determination of any one coupling, 
the values of all the others were held at their Standard Model values. The limits shown 
include both statistical and systematic effects: a) Coefficients of lowest dimension oper- 
ators contributing either to ZZ and ZY production or to Z 7 , ZY and ZZ production; 
b) Coefficients of lowest dimension operators affecting only the V°ZY vertices; c) Coeffi- 
cients of SU(2) x £/(l)-conserving operators according to i) the Gounaris-Layssac-Renard 
constraints and ii) the Alcaraz constraints (see text, section 1.1). 
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of the background level. In the qqvv channel, as in qqqq, the main source of systematic 
error arises from modelling of the qq(j) background, in this case corresponding to the 
kinematic region with large missing energy, and hence low visible qq energy. A study of 
the energy flow in this region using events at the Z peak allowed a determination of the 
effect of this uncertainty in the present analysis; it gives rise to systematic errors in the 
coupling parameters of order 5% - 10% of the values of the statistical errors. Another, 
comparable source of systematic error in this channel comes from the uncertainties in 
the cross-sections for the dominant background channels, particularly Weu production. 
Systematic effects in the l + l^vv channels were found to be negligible. In addition, the 
effects of uncertainties of ±2% in the overall ZZ cross-section and of ±1% in the luminosity 
measurement were considered. The combined effect of all the systematic uncertainties in 
the channels contributing to ZZ production is small, typically ~ 15% of the statistical 
errors, and, as in the case of the /i^-related parameters discussed above, this ratio of 

yOyOyO 

systematic to statistical effects is also applicable to the results for the £ { 1 2 3 m z related 
to on-shell fY parameters shown in table 6a). 

The systematic uncertainties in the study of the qqe + e~ and qqfi + f^~ channels con- 
tributing to ZY production have been described in [9]. Several effects, including uncer- 
tainties in lepton identification, the effect of limited simulated data and, in the qqe + e~ 
channel, identification of fake electrons coming from background channels, combine to 
give a systematic error on the efficiency to select qqe + e~ and qq[i + fi~ events of ±5% and 
a relative uncertainty in the background level of ±15%. In addition, a systematic error of 
±1% in the luminosity measurement was assumed. The overall effect of these systematic 
uncertainties in the determination of the coupling parameters is small in comparison with 
the statistical errors. In combination with ZZ data to produce the results for parameters 
£ 2 ZrT and i^ Zl listed in table 6b), they amount to ~15% and ~5% of the statistical 
errors, respectively. 

In the combination of data from the different final states, Zy, ZZ and Zj*, all the 
systematic effects listed above were treated as uncorrelated except those arising from the 
uncertainty in the luminosity measurement. 

3.2 Discussion 

A few comments may be made on the results shown in table 6 and figures 6-10. 

All the results are compatible with the Standard Model expectation of the absence of 
neutral triple-gauge-boson couplings. The results shown in table 6a) and figures 6 and 7 
demonstrate this conclusion in the effective Lagrangian model of reference [2] for the 
d — 6 operators which, in the on-shell limit, contribute either to ZZ or to Z7 production. 
As mentioned in sections 1.1 and 3.1 (and predicted from studies of simulated events [6]), 
the contribution to these results of the off-shell data included in their determination is 
small: using only the off-shell data leads to precisions poorer by factors of ~ 3 — 7 than 
using the on-shell Z7 or ZZ data. (This effect is observed most strongly in the case of the 
determination of hj, where the interference in the squared matrix element between the 
anomalous and Standard Model amplitudes leads to a relatively precise determination 
of this parameter). Thus these results, with negligible changes, may also be interpreted 
in terms of the parameters hf and jj of on-shell Z7 and ZZ production, listed in the 
right-hand column of the table, and they may be compared directly with other published 
results for these on-shell parameters. 

The results shown in table 6b) and figure 8 examine the possibility of four-fermion 
production via an anomalous V°Z'j* vertex by determining the coefficients of the lowest 
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dimension (d = 8) operators in the model of reference [2] which would contribute to such 
a process. As noted in section 1.1, contributions from these operators affect both the 
Z'-f* and ZZ final states; in the determination of £ 2 Z ' y ' y , the experimental samples from 
the two final states contribute roughly equally to the log likelihood distribution in the 
combination of data, while in the determination of I ZZl the Zq* contribution dominates. 
The results of the fits show that there is no evidence in the data for a CP-conserving 
anomalous coupling at the 7Z7* vertex or for a CP-violating coupling at the ZZ^f* 
vertex. 

The results shown in table 6c) and figures 9 and 10 indicate that there is no evidence in 
the data for SU(2) x £/(l)-conserving anomalous couplings in the models of references [2] 
and [6]. Here again, in the combinations of data from different final states, the contri- 
butions from Z7 production dominate, as can be seen by comparison of the likelihood 
curves of figure 7 and either figure 9 or figure 10, and from the confidence limits shown 
in the table. This arises both because of the sensitivity to hj noted above and because 
of the greater statistical contribution from Z'-f compared to that from ZZ production at 
LEP2 energies. 

4 Conclusions 

A study has been performed of the neutral triple-gauge-boson vertex using DELPHI 
data from the final states Zj, ZZ and Zq* produced at LEP2. The results have been 
interpreted in terms of various models of the interaction Lagrangian proposed in the 
literature. We find no evidence for the production of these states by processes involving 
neutral triple-gauge-boson vertices with either one or two off-shell bosons, nor when the 
data are analyzed in terms of models in which the neutral triple-gauge-boson vertex is 
constrained to be 577(2) x £/(l)-conserving. These conclusions are in agreement with the 
predictions of the Standard Model. 
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Figure 2: a) Distribution of x 7 = i? 7 /i?fc eam , the energy of identified photons normal- 
ized to the beam energy in the data in the uuj channel, summed over all energy points. 
The distribution is shown before imposing the experimental cut at E^ = 50 GeV. The 
experimental data points are shown by dots and the shaded histogram shows the predic- 
tions of the Standard Model for signal plus background. (The background contribution 
is very small, and is not shown separately), b) Distribution of | cosa*|, where a* is the 
decay angle of the quark (or antiquark) in the Z centre-of-mass frame with respect to 
the direction of the Z in the overall centre of mass, for data selected in the q<pi channel. 
The experimental data points are shown by dots, the shaded histogram shows the pre- 
dictions of the Standard Model for signal and background, and the outlined histograms 
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qqVV 

llvv 

background 



Figure 3: Distribution of | cos 9z \ , where 6z is the Z production polar angle, for data 
selected in the ZZ channels. The experimental data points are shown by dots, the shaded 
histograms show the predictions of the Standard Model for the signal and background 
components indicated in the legend, and the outlined histograms the expectations for 
values of /rf = ±1.5. 
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Figure 4: For the qql + l final state: predicted Standard Model distributions of events 
a) in the (M hadrons , M^^-) plane, and b) in the (M hadrons , M e+e -) plane, showing the 
bins used in the fits to the coupling parameters. The sum of all the bins defines the 
Z'-f* sample, c) Expected distribution in the (M q q, M^+ M -) plane, and d) in the (M qq -, 
M e + e -) plane, of the difference between the predictions of the Standard Model plus an 
anomalous contribution, i^ Zl m\ = 3.4, and the Standard Model only. (The parameter 
i ZZl is defined in table 1). Plots a) and b) and, separately, c) and d) were computed 
with the same assumed luminosities. 
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Figure 5: a) Distribution of Mi+i- (I = e,/i), and b) of M q g, for data selected in the 
qqH + H~ and qqe + e~ channels. The experimental data points are shown by dots, the full 
histograms show the predictions of the Standard Model for signal and background, and 
the dotted histograms the expectations when an anomalous contribution, I ZZl m A z = 3.4, 
is present. 



24 




Figure 6: Likelihood distributions for neutral gauge coupling parameters corresponding 
to Lagrangian operators influencing ZZ and Zj* production. The parameters are defined 
in section 1.1; the corresponding on-shell parameters are shown in parentheses on the 
abscissa labels. The distributions include the contributions from both statistical and 
systematic effects. 
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Figure 7: Likelihood distributions for neutral gauge coupling parameters corresponding 
to Lagrangian operators influencing Z7, Z7* and ZZ production. The parameters are 
defined in section 1.1; the corresponding on-shell parameters are shown in parentheses 
on the abscissa labels. The distributions include the contributions from both statistical 
and systematic effects. 
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DELPHI (ZZ,Zy*) 
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Figure 8: Likelihood distributions for neutral gauge coupling parameters corresponding 
to Lagrangian operators affecting only the V°Z / ~f* vertices. The parameters are defined 
in section 1.1. The distributions include the contributions from both statistical and 
systematic effects. 



27 



DELPHI (Zy, Zy*,ZZ) G-L-R SU(2)xU(l) invariance 
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Figure 9: Likelihood distributions for neutral gauge coupling parameters corresponding to 
577(2) x £/(l)-conserving Lagrangian operators satisfying the Gounaris-Layssac-Renard 
(G-L-R) constraints. The parameters are defined in section 1.1. The distributions include 
the contributions from both statistical and systematic effects. 



28 



DELPHI (Zy,Zy*,ZZ) 




-0.2 0.2 



- cot W m z 2 v 2 /4 l g B ( hj Z ) 



Alcaraz SU(2)xU(l) invariance 




-0.05 0.05 

- cot W m z 2 v 2 /4 1 8 A ( h 3 T ) 




-0.2 0.2 



- cot W m z 2 v 2 /4 1 8 B ( h 3 Z ) 



Figure 10: Likelihood distributions for neutral gauge coupling parameters corresponding 
to 577(2) x [/(l)-conserving Lagrangian operators satisfying the Alcaraz constraints. The 
parameters are defined in section 1.1. The distributions include the contributions from 
both statistical and systematic effects. 



